Exact solution and magnetic properties of an anisotropic spin ladder 
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We study an integrable two-leg spin-1/2 ladder with an XYZ-type rung interaction. Exact 
rung states and rung energies are obtained for the anisotropic rung coupling in the presence of 
a magnetic field. Magnetic properties are analyzed at both zero and finite temperatures via the 
thermodynamic Bethe ansatz and the high-temperature expansion. According to different couplings 
in the anisotropic rung interaction, there are two cases in which a gap opens, with the ground 
state involving one or two components in the absence of a magnetic field. We obtain the analytic 
expressions of all critical fields for the field-induced quantum phase transitions (QPT). Anisotropic 
rung interaction leads to such effects as separated magnetizations and susceptibilities in different 
directions, lowered inflection points and remnant weak variation of the magnetization after the last 
QPT. 



I. INTRODUCTION 

Recently the quasi-one-dimensional spin ladder has 
attracted much interest both experimentally and the- 
oretically [1]. More and more ladder-structure 
compounds have been realized, such as SrCu203 
[2], Cu2(C5Hi2N2)2Cl4 [3], (5IAP)2CuBr4-2H20 [4], 
(C5Hi2N)2CuBr4 [5], and so forth. Although many lad- 
der compounds can be well described by simple isotropic 
ladders, the structural distortion and the spin-orbital 
interaction of the transition ions can lead to various 
magnetic anisotropies. Besides the spin-orbital inter- 
action, both on-site Coulomb exchange interaction [6,7] 
and nonlocal Coulomb interaction [8] also influence the 
anisotropy. Anisotropic interaction from bond buckling 
has been recently found in copper-oxide ladder com- 
pounds CaCu203 [9] due to an angle deviation from 180° 
in the Cu-O-Cu bond [9-ff]. An anisotropic rung inter- 
action was considered in Ref. [12] motivated by CaCu203 
[9] and a two-leg spin ladder with an XXZ-rung interac- 
tion was derived in the presence of the Dzyaloshinskii- 
Moriya interaction and the Kaplan-Shekhtman-Entin- 
Wohlman-Aharony interactions. When the Cu-O-Cu 
bond is near 90°, the rung interaction is weak in the 
copper-oxide ladder. Spin anisotropy in the exchange 
interaction also exists in strongly-coupled ladder com- 
pounds such as (pipdH)2CuBr4 [f3]. On the other hand, 
real spin ladder compounds are usually described by the 
standard Heisenberg ladder model, which is not exactly 
soluble, turning the computation of the physical proper- 
ties for the ground state (GS), the gap, the thcrmody- 
namical quantities and other relevant properties in the 
presence of temperature and magnetic fields, rather dif- 
ficult. Usually, just numerical calculations and pertur- 
bative schemes can be applied. Recently it was shown 
that the integrable spin ladder model [f4] can describe 
the properties of strongly-coupled spin ladder compounds 
when a rescaling parameter is introduced [15, f6]. There- 
fore, it can be expected that integrable ladders with 



anisotropic rung interactions can provide some meaning- 
ful information in the physics of anisotropies. In addi- 
tion, anisotropic rung interactions also provide us with 
some more adjustable parameters that may be useful in 
fitting the experimental data of compounds. 

In the present paper we shall consider the anisotropy 
in the rung interaction and the corresponding magnetic 
anisotropic effect by solving an integrable spin ladder 
with a general XYZ rung interaction. By means of the 
thermodynamical Bethe ansatz (TBA) and the high tem- 
perature expansion (HTE) [15,17,18], we investigate the 
influence of the anisotropic rung interaction on the quan- 
tum phase transitions (QPT) and the magnetic proper- 
ties. The contents are arranged as follows: (i) In the 
second section we present the model and the exact rung- 
state basis in the presence of a magnetic field. Then 
the model is solved by the Bethe ansatz (BA) approach, 
(ii) The third section gives the TBA equations for the 
GS and the HTE for physical properties at finite tem- 
peratures, (iii) In the forth section, we study the field- 
induced QPT's. The rung anisotropy provides two kinds 
of gapped ladders, respectively with one and two compo- 
nents in the GS. The analytic expressions are obtained 
for all the critical fields of the corresponding QPT's. The 
rung anisotropy also leads to a separation of the magne- 
tizations and susceptibilities in different directions. The 
magnetization inflection point (IP) may be lowered from 
the half-saturation and in the two-component gapped 
ladder the IP is even not invariant under different tem- 
peratures. A remnant variation of magnetization can be 
found after the last QPT. In the fifth section we give a 
summary of our results. 



II. THE MODEL, EXACT RUNG STATES AND 
BA SOLUTION 

We shall consider a spin-1/2 two-leg spin ladder model 
with a general XYZ-type anisotropy in the rung interac- 
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tion, whose Hamiltonian reads 

n = no + nxYz + M, 

Jo ' 



i=l 



i 

M = -gHj2iSt+Tn^ (1) 



where S and f are the spin operators for the two legs 

and g is the Landc g factor in the direction of the field. 
Jo is the average of the leg interaction and 7 is a rcscal- 
ing parameter. The value j — 4 was introduced in Rcfs. 
[16,15] in fitting with isotropic ladder compounds in the 
presence of weak rung. The bulk part Ho with the permu- 
tation operator P,,,+i = (25^, • + i)(2f, • f,+i + i) 
exhibits the SU(4) symmetry [19]. Isotropic integrable 
spin ladder [14] has identical Jx, Jy and Jz- 

When the rung interaction is strong, it is favorable 
for the spin ladder system to form rung states, since the 
leg interaction is too weak to take apart the rung state. 
Anisotropy in the rung interaction leads to the collapse of 
the conventional singlet and triplet rung states from the 
isotropic ladder, even in the absence of the field. How- 
ever, we find a new exact basis valid both in the absence 
and presence of an external magnetic field. 
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where 



Then the corresponding rung energies include the Zee- 
man energy in a nonlinear way, 

El = —^{Jx + Jy + Jz), 
E2 = -^{Jx + Jy - Jz), 



Ei 



{9H)' + j^iJx-Jy)' + ^Jz. 



(4) 



The rung states {(pi \ i = l,---,4} provides a new 
fundamental representation of the SU(4) Lie algebra 
Smfi = ^n,i'Pm with Commutation relations of the gen- 
erators [5'"j,S'[] = (5„_fe5'„ — Sjn^iSJ^. Based on this 
SU(4) realization and the vanishing commutation rela- 
tions [HxYZ,'Ho] = [AdjHo] = 0, one can solve the 



model (1) via the BA approach [20]. The BA equations 
arc the same as those obtained for the SU(4) model [21] 
and for the SU(3)®U(1) spin ladder [14]. Here we present 
the BA equations together with the eigenenergy 
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E = -Y^ 2nai{iJif^) + EiNi, (5) 
i=i 



i=l 



where S^(/x^''„ 
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0, and 1 < A: < 3; 
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. There are L rungs, Ni is the total 



number of rung state (fii and /z^- is the rapidity. M^''^ is 
the total rapidity number in the fc'th branch. 

From the BA equations and the eigenenergy we can 
apply the TEA and HTE to study the collective proper- 
ties of the model. When the field is applied in the x or y 
direction, one only needs to permute the anisotropy pa- 
rameter values {Jx, Jy, Jz} as well as the corresponding 
g factors. We incorporate g into the field unit in plot- 
ting figures to investigate the net effect of the anisotropic 
rung interaction. For a powder sample, a simple way is 
to take an average of the three directions. 



III. TBA AND HTE 

By adopting the string conjectures [22] and applying 

the Yang- Yang method [23] at the thermodynamic limit, 
one can obtain the GS equations for three dressed ener- 



gies e 



), (6) 



where e^°^ = e^^^ = and the symbol * denotes the con- 
volution. The basis order is chosen as (ipp^cpp^ipp^ipp^)'^ , 
where Pi <E {1,2,3,4} and ipp-^ is energetically the most 
favorable state while (pp^ is the least favorable one. For 
the chosen order the driving term is given by 17^'^ = 
i?p.^j — Ep.. The GS is composed of Fermi seas filled 
by negative e^*^". If some branch of the dressed energy 
is all positive, then the corresponding excitations to this 
branch is gapful. A QPT occurs at the point where the 
gap is closed. We shall apply these TBA equations to 
analyze the field-induced QTP for the GS. 

For the finite temperature case the TBA involves an 
infinite number of coupled integral equations. In the 
present paper we shall apply the HTE [18,17,15] from 
T-system [24] within the Quantum Transfer Matrix for- 
malism [25] , which involves only a finite number of inte- 
gral equations and consequently is more convenient. Fol- 
lowing Refs. [15,18,17], one can obtain the free energy 
/ for per rung in high temperatures. Here we present 
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the first four terms which dominate the physics for high 
temperatures: 



/ = -T ( Co + C,{^) + + C3(4^' 



(7) 



where T is the temperature, the coefficients are 

_ lOQ 18Q^ 80Q^ 8Q- 24QQ_ 4 

with the definitions 

Q = 2cosh(i/3J^) + 4cosh(i/?Ja:+j,) cosh(/3/i), 
g± = 2e(±'5''^/4) cosh(i^J, +y) + 2e'^^'^'^''^'^^ cosh(/3/i), 

/i = v^(p7TXI77l6, J.±y = J. ± and /? = 1/T. 
One can get higher orders for lower temperatures. The 
magnetization and the susceptibihty can be easily ob- 
tained by M = -df/dH, x = OM/dH. The rescahng 
parameter 7 = 4 in the isotropic case [16,15] gives the 
leading term of the gap for the integrable ladder in fitting 
with the experimental ladder compounds. If the rung Jo 
is weak, then the HTE gives a valid result even for low 
temperatures due to the large rescaling 7. 



IV. PHASE TRANSITIONS AND MAGNETIC 
PROPERTIES 

A. One-component gapped ladder 

For different anisotropics there are two different 
gapped ladders. In one case, only one component ipi 
exists in the gapped GS when H = 0. In the other case 
both (fi and ip2 are involved in the gapped GS. First we 
discuss the former case which happens more likely. It 
requires 



Jz + Ja- \ Jx - Jy\ > 16Jo/7> 

+ Jy> 8J0/7, 



(8) 
(9) 



where J3 = Jz + Jy + Jz- Condition (8) expels the 
components ip^ and (^4 from the GS, while condition 
(9) excludes the component ip2. The field will bring 
(fis down to the ground state and close the gap A = 
min{£J2, E3, £4}— i^i— 4 J0/7 at a critical field Hd, which 
leads to the first quantum phase transition (QPT). Fur- 
ther increase of the field will bring all components of ipi 
out of the GS and another gap A = Ei — E3 — 4J0/7 
opens at the critical field Hc2, which characterizes an- 
other QPT. The factor 4J0/7 in the gap comes from the 
maximum depth of the first dressed energy branch. It is 
easy to see that only the components ifi and ^3 compete 



in the GS involving one branch of dressed energy, since 
the GS only consists of ipi in the absence of the field 
while only (p^, is lowered in energy level when the field 
is applied. The analytic expressions of two critical fields 
can be obtained exactly as 



Hcl = Jzh + J.Jy + 64(^)2 - 8^ (J, + J3), 

Hc2 = Ygsj J-J^ + J-^y + 64( + 8 Y (-^^ + ■^3)- (10) 

Setting = Jy = J_j, recovers the result for isotropic case 
[16], as expected. A weak anisotropy will lead to differ- 
ent critical fields and consequently separate the magne- 
tizations in different directions. We give an example of 
the magnetization with weak anisotropy in Fig.l, a low- 
temperature magnetization was presented for comparison 
in the z direction. The corresponding low-temperature 
magnetizations for all three directions are presented in 
Fig. 2, obtained from the HTE. Magnetizations in differ- 
ent directions for strong anisotropy are demonstrated as 
an example in Fig.3 for the GS and in Fig. 4 for a low 
temperature. 

Before the gap is closed at Hd, the gap A near Hd 
can be expanded to a simpler form 



A ^ ci {Hd - H) , 



(11) 



ci = g^Hdl \j {gHdY + ^(Jcc - JyY- Considerable ex- 
citations can be stimulated by the temperature T if T 
is in the order of the gap T ^ (Hd — H), the magne- 
tization will rise from zero before the field reaches the 
critical point. An expansion based on small Fermi points 
[16] gives the zero-temperature critical behavior in the 
vicinity of Hd 
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1/2 



(12) 



Here {M^)^ is the magnetization of a single rung state 
(fs, it also varies with the field due to the anisotropic 
rung interaction, as we will discuss below in (17). For 
the lowest order in the critical behavior, {M'')^ takes the 
value at the critical point Hd- This oc (if - HdY^"^ 
critical behavior, typical for gapped integer spin antifer- 
romagnetic chains [26], is buried by the afore- mentioned 
temperature effect. This temperature effect can be seen 
in Fig. 5, in which the magnetization of the z direction 
at T = O.5J0 becomes considerable at the field H = 
Hd — O.SJq. Actually the magnetization at T = 0.5 Jo 
increases nearly in a linear way before the Hd- 

A special point in the magnetization is the inflection 
point (IP) Hip, which is an invariant point under low 
temperatures, 



gH. 



IP 



\yJ{Jz + J.){Jz + Jy) (13) 
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where the two components (pi and 953 have the same rung 
energies Ei = and the same proportion A'^i = A^3 in 
the GS. The excitations to Lp2 and (/?4 are gapful, the gap 
can be obtained exactly from (6) 



A/p = min{A/P2, A/P4}, 
A/P2 = (Jx + Jy)/2 - (21n2)Jo/7, 
A/P4 = (J^ + J3)/2 - (2 In 2) Jo/7- 



(14) 



At low temperatures, the excitations to (^2 or ^4 are dif- 
ficult to stimulate, while the temperature does not influ- 
ence the relative proportion between ipi and 993 due to 
their same energies at the IP. Consequently the propor- 
tions of ip\ and ifz remain almost unchanged when the 
temperature varies. Therefore the magnetization at Hjp 
also remains unmoved when the temperatures changes, 
and the magnetization curves of various temperatures 
cross each other at the same point M/p, which is shown 
by the curves for temperatures T = 0, 0.5Jo, 0.75Jo and 
Jo in Fig. 5. This requires low temperatures 



T«; Ajp 



(15) 



as well as the gapped ladder conditions (8) and (9), 
A/p is the excitation gap to 1^32 or ^p^ in (14). When 
the temperature is sufficiently high such that the exci- 
tations to (^2 or ipi arc considerable, the involvement of 
these components reduces the proportion of (^3 which has 
the highest magnetization. The components (^2 and (^4 
have zero and negative magnetizations, respectively. As 
a result, the magnetization at iJ/p deviates from M/p 
and move downwards. We show this moving by mag- 
netization curves at temperatures T = 10 Jo, 20 Jo in 
Fig. 5, for which one can find observation examples in 

CU2(C5Hi2N2)2Cl4 [3]. 

The magnetization at the IP can be worked out as 
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IP 



HjpH 
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(+)' 



(16) 
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where H^jp = ^Hjp + (J, - JyY/iAgY. 



For the 

isotropic ladder, h'"^ = Hjp and consequently Mfp 
is located at the half of the saturation magnetization 
= g [15]. The anisotropy lowers the magnetization 
of the IP due to ijj+^ > Hip, i.e. Mfp /Ml < 1/2. 
Physically, the anisotropy in x,y directions hybridizes the 
elemental state |||) into (^3 so that (f^ is not a pure fully- 
polarized elemental state |tT) as in the isotropic case. For 
XXZ-type rung interaction Mjp is half-saturation when 
the field is oriented in z direction, but also lowered when 
the field is applied in other directions. This IP lowering 
effect is more obvious for the strong anisotropic case, we 
give an example in Fig. 3. As one can see from this fig- 
ure, besides the strong separation of the magnetization 
in different directions, the IP points in x and z directions 
move below the half saturation point. 

In addition to the separation of the magnetization in 
different directions and the lowering of the inflection 



points, another property in the anisotropic case is the 
remnant variation of the magnetization after the second 
phase transition. The fact that magnetization increases 
between Hd and iJc2 mainly comes from the propor- 
tional competition between the two state (pi and (^3, i.e., 
more rungs are occupied by Lps when the field gets higher. 
The single-rung magnetization in state (^3 can be ob- 
tained explicitly 



{M% = g 



1 



(17) 



where r] increases with the field from the expression (3). 
The long-dashed line in Fig. 3 gives an example of (M^)^ 
in z direction, which increases from zero from the be- 
ginning of the application of the magnetic field. If Hd 
is small, then the increment of {M^)^ also contributes 
with an important part to the growth of the magneti- 
zation. Otherwise for higher Hci, the change of (M^)^ 
contributes less to the growth of the total magnetiza- 
tion, since {M^)^ has slowed down in increasing before 
the first quantum phase transition occurs. However, the 
competition between (fi and (^3 comes to end after the 
second phase transition and the magnetization is com- 
pletely (M^)^. This gives a remnant variation of mag- 
netization even after the second phase transition, since 
{M^)^ is still approaching to the saturation limit. This 
remnant variation of magnetization is illustrated for the 
ground state in Fig. 3 and can also be seen for the tem- 
perature case ( Fig. 4 ). 

Examples of the magnetic susceptibility in the three 
directions are plotted in Fig. 6 for weakly anisotropic 
rung and in Fig. 7 for strongly anisotropic rung. Weak 
anisotropy separates the heights of the magnetic suscep- 
tibility peak, while a strong anisotropy also leads to an 
obvious separation of the whole susceptibility including 
the peak positions. 



B. Two-component gapped ladder 

Anisotropy in the rung interaction also provides an- 
other possibility of a gapped ladder, in which not only the 
rung state ipi but also ip2 are involved in the GS before 
the field is applied and brings about the first QPT. The 
single-state energy difference is E2 — Ei = {J^ + Jy)/'2- 
The larger is the difference, the more strongly ipi and 
(f2 will expel each other in the Fermi sea. The two- 
component gapped ladder requires 



\Jx + Jy\< 8J0/7, 



(18) 



so that (fi and (f2 are close enough in the energy levels 
to exist in the GS at the same time in the absence of the 
field. Also a strong J^ is needed to expel cp^ and (f4 from 
the gapped GS before the field is applied, approximately 

J, > 41n2^ + i I J, - Jy\ + —L^ (J, + Jyf . (19) 
7 z on Jo 
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For simplicity, we assume + Jy > so that tpi has 
lower energy than <p2, one only needs to change + Jy 
to be —{Jx + Jy) for lower (^2- The first QPT occurs when 
the field brings down (^3 and gets involved in the GS, the 
critical field can be obtained with the help of the Wiener- 
Hopf technique [27] which is valid for large Fermi points 
(Fermi surface in one dimension). Explicitly we have 



cl 



2 In 



167r2 J0/7 
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(20) 



which gives a good approximation if the value of Jx + Jy 
is not very close to 8 J0/7. Further increase will lower the 
energy of (^93 below ipi and and bring them out of the 
GS one by one. The component variations in the QPT 
are {iyJi,<p2} ^ ¥'2, (ps} ^ ^ W^}, where 

each arrow indicates the occurrence of a QPT. Since (^2 
has also zero magnetization, the total magnetization also 
remains null in the gapped phase before the first QPT. 
The zero-magnetization component <^2 gets out of the GS 
after Hjp if 



while for 



Jx + Jy< (4 In 2) Jo /7, 



(4 In 2) Jo /7 < Jx + Jy < 8J0/7, 



(21) 



(22) 



(p2 is brought off the GS before Hjp. These happen at 
the second QPT with an approximate critical field 



gH, 



c2 



^x—y 

(23) 



where 5 = Jx+y — 41n2(Jo/7). The expression (23) can 
give a satisfactory approximation when the value oi \5\ 
is not near 4J0/7. The exact critical field Hc3 for the 
third QPT is the same as Hc2 in (10). When the ex- 
ample in fig. (8) has numerical points Hd — I.O6IJ0 and 
iJc2 = I.232J0, the expressions (20) and (23) provide 
analytic results Hd = 1.064Jo and iJc2 = 1.231 Jo. 

The IP in the one-component gapped ladder case will 
not be invariant in the two-component ladder case. If 
the component <^2 gets out of the GS after Hjp, the (^2 
is gapless. Although the components (pi and Lp:^ still have 
the same proportion at the IP, any small temperature will 
excite more components of f2 and consequently decreases 
the proportion of (pi and (p^. Therefore the temperature 
will lower the total magnetization from that of the GS. 
If the component ip2 gets out of the GS before iJ/p, with 
the condition (22), the IP is hardly an invariant. Despite 
of the existence of a gap for excitations to <^2 at Hjp, 
the gap is actually quite small 



relative to the strong rung interaction. So a low tem- 
perature of order Jq will still stimulate considerable ex- 
citations to (p2 and lower the magnetization at i?/p. We 
illustrate this by an example in Fig. (8). 



V. SUMMARY 

We have introduced a two-leg spin-1/2 ladder with a 

general anisotropic XYZ rung interaction. In particular, 
the exact rung state basis for this model was found. We 
have studied the effect of the anisotropic rung interaction 
by solving the integrable ladder in the context of the ther- 
modynamical Bethe ansatz and the high-temperature ex- 
pansion. Two kinds of gapped ladders were provided, re- 
spectively involving one and two components in the GS in 
the absence of the magnetic field. We have obtained ana- 
lytically all the corresponding critical fields for the field- 
induced quantum phase transitions. The magnetizations 
and susceptibilities in difi'erent directions separates under 
the rung anisotropy. The magnetization inflection point 
is lowered from the half-saturation and a weak chang- 
ing in magnetization still remains after the last quan- 
tum phase transition. The inflection point in the two- 
component gapped ladder case is not invariant as in the 
one-component gapped ladder case due to field-deduced 
three-component competition or small excitation gap. 
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FIG. 1. Magnetization versus the magnetic field at zero temperature for the one-component gapped ladder with a 
weakly-anisotropic rung interaction, = 8Jo, Jy = 9Jo, Jz = 10.5Jo, with the rescaling 7 = 4. Ms is the saturation 
magnetization. In the gapped phase H < Hci, only one component ipi exists in the ground state. To study the net effect 
of the anisotropic rung interaction, we incorporate the g factor into the field. The weak anisotropy in the rung separates the 
magnetization in different directions. The zero temperature magnetization is obtained from the thermodynamical Bethe ansatz 
(TBA). Also for comparison with the finite temperature case, a magnetization at T = 0.5 Jo obtained from the high-temperature 
expansion (HTE) is presented in the z direction. 
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FIG. 2. Magnetization versus the magnetic field for T = 0.5Jo for tiie one-component gapped ladder with the weak anisotropy, 
Jx ~ 8 Jo, Jy — 9 Jo, Jz ~ 10. 5 Jo and "f — 4. The magnetizations are obtained from the HTE, which coincides with the 
magnetization separation in the zero temperature case obtained from the TBA. 
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FIG. 3. Magnetization versus the magnetic field at zero temperature for the one-component gapped ladder with a 
strongly-anisotropic rung interaction, — 3Jo, Jy = 15Jo, Jz = 6Jo, 7 = 4. The strong anisotropy leads to a strong 
separation of the magnetization. The inflection point (IP) is lowered from the half saturation. A temperature magnetization 
from HTE is present to demonstrate the IP. Note that the magnetizations do not reach the saturation after the second quantum 
phase transition at Hc2, there still remains a weak variation of magnetizations. This remnant magnetization variation comes 
from the single-state magnetization M| of ^3 which is a mixture of full-polarized states ||t) and the lowest-magnetized state 
III). The variation of Aff is illustrated by the long-dashed line for the whole process of the field application. 
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FIG. 4. Magnetizations versus the magnetic field for different directions in temperature case for tlie one-component gapped 
ladder with the strong anisotropy, — 3 Jo, Jy — 15 Jo, Jz = 6 Jo, 7 = 4. 
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FIG. 5. Magnetization versus the magnetic field at different temperatures for the strongly-anisotropic rung, = 3 Jo, 
Jy = 15Jo, Jz = 6Jo, 7 = 4. Low temperature magnetizations in T = 0, 0.5Jo, 0.75Jo and Jo cross the inflection point (IP). 
Higher-temperature magnetizations at T = 5Jo and lOJo do not go through the IP, as the gap for excitation to <^2 is overcome 
by the temperature stimulation. 
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Magnetic susceptibility against the temperature in different directions for the weakly-anisotropic case, J^ = 8 Jo, 
Jz = 10.5Jo and j — 4. The weak anisotropy separate the heights of the susceptibility peaks. To see the net effect 



of the rung anisotropy, we plot the figures using the same g factors for the three directions. 
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FIG. 7. Magnetic susceptibility versus the temperature in the strongly-anisotropic case, Ja; = 3Jo, Jy = 15 Jo, Jz — 6 Jo, 
J — 4. The strong anisotropy separate not only the peak heights but also the whole shape including the peak positions. 
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FIG. 8. Magnetizations in z direction of a two-component gapped ladder case, — Jy — 0.5 Jo, Jz = 5 Jo, 7 = 4. The 
curves a, b, c, d and e are plotted in different temperatures, respectively, T= 0, 0.5Jo, Jo, 1.5Jo and 5Jo. In the gapped phase 
H < iJci, two components ipi and ip2 involve in the ground state (GS) (T=0). The component ips begins to enter the GS at the 
critical point Hd and reaches the same energy as </P2 at H23- The component ip2 get out of the GS at the second critical field 
Hc2- The magnetization curves of temperatures do not go through the IP as in the one-component gapped ladder case. When 
the field is applied in a; or y direction the GS magnetization will increase from the beginning due to the gapless excitation in 
these directions. 
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